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We introduce a class of models defined on ladders with a diagonal structure generated by rip 
plaquettes. The case Up — 1 corresponds to the necklace ladder and has remarkable properties 
which are studied using DMRG and recurrent variational ansatzes. The AF Heisenberg model 
on this ladder is equivalent to the alternating spin-l/spin-i AFH chain which is known to have 
a ferrimagnetic ground state (GS). For doping 1/3 the GS is a fully doped (1,1) stripe with the 
holes located mostly along the principal diagonal while the minor diagonals are occupied by spin 
singlets. This state can be seen as a Mott insulator of localized Cooper pairs on the plaquettes. A 
physical picture of our results is provided by a tp — Jp — td model of plaquettes coupled diagonally 
with a hopping parameter t^. In the limit td —> oo we recover the original t — J model on the 
necklace ladder while for weak hopping parameter the model is easily solvable. The GS in the 
strong hopping regime is essentially an "on link" Gutzwiller projection of the weak hopping GS. 
We generalize the tp — Jp — td model to diagonal ladders with rip > 1 and the 2D square lattice. 
We use in our construction concepts familiar in Statistical Mechanics such as medial graphs and 
Bratelli diagrams. 

PACS number: 74.20.Mn, 71.10.Fd, 71.10.Pm 



I. INTRODUCTION 

Ladders provide a class of interesting theoretical mod- 
els for studying the behavior of strongly correlated elec- 
tron systems. Besides representing simplified models for 
actual materials, ladders offer a possible way of interpo- 
lating between 1 and 2 spatial dimensions with the hope 
that they will yield insights into the physics of 2D sys- 
tems, such as the Cu02 planes of the cuprates (for a 
review see 

It has been found that ladders exhibit quite different 
behavior depending on whether the number of legs ng is 
even or odd. Antiferromagnetic spin ladders with ng odd 
are gapless with spin-spin correlation functions decaying 
algebraically, while even leg ladders are gapped with a 
finite spin correlation length. Upon doping, these two 
types of ladders also behave differently concerning the 
existence of pairing of holes or spin-charge separation. 
In the limit where the number of legs goes to infinity 
the spin gap of the even spin ladders vanishes exponen- 
tially fast, in agreement with the gapless nature of the 
2D magnons [§. On the other hand, the antiferromag- 
netic long range order (AFLRO) characteristic of the 2D 
antiferromagnetic Heisenberg (AFH) model can be more 
naturally attributed to the quasi long range order of the 
odd leg ladders. It thus seems that one has to combine 
different properties of the even and odd, doped and un- 
doped ladders in order to arrive at a consistent picture of 



the 2D cuprates. Ladder systems are sufficiently interest- 
ing on their own to deserve detailed studies, in addition 
there are a variety of materials which contain weakly cou- 
pled arrays of ladders |^ . 

In this paper we study of a class of ladders character- 
ized by a diagonal structure which provides an alternative 
to the aforementioned route to 2D. We shall call these ob- 
jects diagonal ladders in order to distinguish them from 
the more familiar rectangular shaped ones. Diagonal lad- 
ders are labelled by an integer Up = 1,2,... which gives 
the number of elementary plaquettes needed to generate 
the entire structure. The first member of this family, 
i.e. Hp — 1, is also known as the necklace ladder and 
it consists of a collection of N plaquettes joined along 
a common diagonal. In this paper we shall focus on 
the necklace ladder, although the other cases will also 
be briefly considered. 

The original motivation of this work was to under- 
stand the fully doped stripes in the (1,1) direction 
that have been observed experimentally in materials like 
Lai_2,Sr2,Ni04 in Hartree-Fock studies of t — J and 
Hubbard models jsj, and numerically in density matrix 
renormalization group § (DMRG) studies of the t - J 
model [Q . The simplest possible toy model of this type 
of stripes is provided by a necklace ladder with a hole 
doping equal to 1/3. As we shall see this doping plays 
an important role in our work. 

Lattices similar to the diagonal ladders, but with ad- 
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ditional one-electron hopping terms along the major and 
minor diagonals of each plaquette, and with J/t = 
0, have been solved exactly for certain fillings 
Giesekus [To| ] has shown that for the corresponding ver- 
sion of the necklace ladder, in the case when all of the 
one-electron hopping terms are equal and the hole dop- 
ing is set to X = 1/3, the model has a short range RVB 
ground state, in which the static correlations exhibit an 
exponential decay and the dynamic correlation functions 
exhibit a gap in their spectral densities. 

Let us also note in passing that diagonal ladders have 
recently appeared as constituent parts of some interesting 
materials like Sro.4Ca13.6Cu24O41.84 known for its super- 
conducting properties at high pressure [pl| 

The organization of the paper is as follows. In Section 
II we define the diagonal ladders from a geometrical view- 
point and compare them with the more familiar ladder 
structures. In Section III we study the AF Heisenberg 
model of the necklace ladder. In Section IV we study the 
t — J model on the necklace ladder and show the con- 
servation of the parity of the plaquettes. In Section V 
we present the ground state (OS) structure of a necklace 
t — J ladder with 7 plaquettes, obtained with the DMRG 
and recurrent variational ansatz (RVA) methods. In Sec- 
tion VI we study in more detail the structure of the GS 
at doping 1/3. In Section VII we introduce a generalized 
t — J model on an enlarged necklace ladder, called the 
tp — Jp — td model, and use it to give a physical picture of 
the results of Sections V and VI. In Section VIII wc define 
the tp — Jp — td model on diagonal ladders with more than 
one plaquette per unit cell and on the 2D square lattice. 
In Section IX we state our conclusions. There are three 
appendices which give the technical details concerning 
the RVA calculations (Appendix A), the complete spec- 
trum of the t — J Hamiltonian on a plaquette (Appendix 
B) and a plaquette derivation of the equivalence between 
the spin 1 AKLT state of a chain and the dimer-RVB 
state of the 2-leg AFH ladder (Appendix C). 




(b) 



(c) 



FIG. 1. Examples of diagonal ladders with a number of 
plaquettes Up = 1, 2, 3 in the unit cell. 



II. GEOMETRY OF DIAGONAL LADDERS 

A diagonal ladder can be characterized by the number 
of plaquettes np of the unit cell and the number N of 
these cells. In Fig. 1 we show diagonal ladders with 
Up — 1,2 and 3. There are Up + 2 sites per unit cell. 
Assuming open boundary conditions the total number of 
sites is then given by Ng = {up + 2)N + np. 



(a) (b) 

FIG. 2. (a)Weak coupling (dashed lines) regime in rect- 
angular ladders, (b) Strong couping (dashed lines) regime in 
rectangular ladders. 




FIG. 3. Plaquette construction (right) of a diagonal ladder 
(left). 

Rectangular ladders can be seen either as a collections 
of legs coupled along the rungs or as collections of rungs 
coupled along the legs (see Fig. 2). This geometric fea- 
ture is the basis of the weak coupling and strong cou- 
pling approaches to the various physical models defined 
on ladders. Thus for example the Heisenberg model on 
the ng-leg ladder is usually defined with an exchange cou- 
pling constant J|| along the legs and an exchange coupling 
constant J± along the rungs. The weak and strong cou- 
pling limits correspond to the cases where J|| >> Jl and 
J|| << J_L respectively. 

On the other hand diagonal ladders do not admit such 
a simple construction. The most natural interpretation 
is to regard them as collections of plaquettes joined along 
their common diagonal (see Fig. 3). The trouble with 
this construction is that it does not preserve the number 
of sites! Indeed one has to fuse the points on the principal 
diagonal of the plaquettes before getting the actual neck- 
lace structure. We shall resolve this problem in Section 
VII on physical grounds. 
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III. THE SPIN NECKLACE LADDER 

Let us begin by considering the AFH model on the 
necklace ladder of Fig. 1(a). The Hamiltonian of the 
model is simply, 



H 



S7" • S^" 



(1) 



<«j> 



where J is a positive exchange coupling constant and the 
sum runs over all links < i, j > of the ladder. We shall 
label the sites of the n}^ plaquette as in Fig. 4. The 
Hamiltonian (lil) then becomes 



N 



H — J (Sl,„ + S2,n) • (Sa^n + S3_„_i) 



(2) 



where Sa,n is a spin 1/2 operator acting at the a = 1, 2, 3 
position of the rfi^ plaquette. Eq.(|^) implies that H de- 
pends on the spins of the minor diagonals through their 
sum 



Si .,1 + S2,n 



(3) 



(l,n) 



(3,n-l) 




(3,n) Plaquette 



(2,n) 



FIG. 4. Diagonal coordinates for a single plaquette. 

At this stage we are free to choose the spins of the n}^ 
diagonal in the singlet (5*12, « = 0) or triplet ( 5i2,,i = 1) 
representations. In the latter case the Hamiltonian (||) 
becomes entirely equivalent to that of an alternating 
spin-l/spin-^ chain. Choosing a singlet on the minor di- 
agonal of a given plaquette amounts to adding a spin zero 
impurity on the corresponding spin-1 site in the alternat- 
ing chain, which therefore breaks into two disconnected 
pieces. The net result is that the spin necklace ladder in 
fact describes alternating spin-l/spin-i chains with all 
possible sizes. 

Fortunately, the alternating spin-l/spin-i chain has 
been the subject of several studies concerning the ground 
state (GS), excitations, thermodynamic and magnetic 
properties p^-p^. The GS turns out to be ferrimag- 
netic with total spin given by sq — N/2 where N is the 
number of unit cells of the chain. There are gapless exci- 
tations to states with spin sq — 1 and gapped excitations 



to states with spin sg + 1- In spite of the existence of gap- 
less excitations, the chain has a finite correlation length 
^ ~ 1, defined from the exponential decay of the spin- 
spin correlation function < Si • Sj > after subtraction 
of the ferrimagnetic long range contribution. These re- 
sults have been obtained by a combination of spin-wave, 
variational and DMRG techniques, with very satisfac- 
tory quantitative and qualitative agreement among them 
PIP- 




t 



0.45 



FIG. 5. DMRG results for the spin configuration for the 
G.S. of an undoped 8x3 necklace ladder. The state has total 
spin So — 4. The length of the arrow is proportional to (Sz), 
according the scale in the box. 

We have confirmed some of these properties by apply- 
ing DMRG and variational methods to the spin neck- 
lace ladder. In Fig. 5 we present a snapshot of the 
spin configurations of the GS of an 8 x 3 ladder, ob- 
tained with the DMRG, which has total spin sg = 4. 
Wc find that the mean value of the spins near the center 
of the system are given by < 5^ „ >=< 5*1 „ >= 0.396 
and < 6*1 >= —0.292 in agreement with the results 
of Ref. |l|j namely < „ >= 0.39624 and < >= 
—0.29248. Also using variational RVA methods we have 
found < Sl„ >= 0.4160 and < S^^ >= -0.2893. 



+ V 




2N-2 > f ^ t 



FIG. 6. Diagrammatic representation of the Recurrent Re- 
lations generating the G.S. of an undoped necklace ladder 
using the variational RVA method (see Appendix A) . 

The existence of a very short correlation length sug- 
gests that the ferrimagnetic GS is an adiabatic deforma- 
tion of the Neel state, which can be described by a short 
range variational state. References |jl^,|l^ propose sev- 
eral variational matrix product states ]15|]. It is more 
convenient for our purposes to use the Recursive Vari- 
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ational Approach (RVA) of references |lJ,|T^, in order 
to deal with doped and undoped cases on equal footing. 
The GS of a ladder of length N is built up from the 
states with lengths — 1, — 2 and eventually A^ — 3 
if A^ is odd. The GS thus generated is a order RVA 
state. In Fig. 6 we show a diagrammatic representation 
of the corresponding recurrence relations (we leave for 
Appendix A the technical details). The GS energy per 
site of the associated alternating chain that we obtain is 
given by — 0.7233J, which is to be compared with the 
extrapolated DMRG results —0.72704 J or the spin-wave 
value -0.718 J of Pati et. al. 111. 



IV. THE t ~ J MODEL ON THE NECKLACE 
LADDER 

The Hamiltonian of the t — J model is given by 





+ h.c.) Vg 



(4) 



where the Ci^s ( c| s ) is the electron destruction (creation) 
operator for site i and spin s, Ui is the occupation number 
operator, and Vg is the Gutzwiller projection operator 
which forbids doubly occupied sites. The density-density 
and kinetic terms in (^) can be written in a form similar 
to (^) for the exchange terms. This suggest that there will 
also be a "decoupling" of degrees of freedom associated 
with the transverse diagonals. The simplest way to see 
how this decoupling works is as follows. 

For the necklace t — J ladder, there is a parity pla- 
quette conservation theorem |l0|] : the t — J Hamiltonian 
on a necklace ladder commutes with every graded per- 
mutation operator P„ associated with the minor diago- 
nal of the n*'' plaquette. Here the permutation operator 
P„{n — 1, . . . , N) is defined by its action on the fcrmionic 
operators, which is trivial at all the sites except for those 
on the minor diagonal of the n^^ plaquette where it acts 
as 



Pn C(l,„),s Pji = C(2,„),s 
Pn C(2,n),s Pn — '^(l,n),s 



(5) 



Of course the spin and the density number operators at 
the sites (1, n) and (2, n) are also interchanged under the 
action of P„. The above theorem is the statement that 
Pn commutes with iJf.j, Eq.(4), for all n, 



[Ht,j, Pn] = for n = l, 



(6) 



and can be easily proved. Eq. (^ is not special to the t—J 
Hamiltonian, since any other lattice Hamiltonian having 
the permutation symmetry between the two sites on the 
minor diagonal of every plaquette would share this same 
property. 

The immediate consequence of (^) is that we can si- 
multaneously diagonalize the Hamiltonian Ht^j and the 
whole collection of permutations operators P„, whose 
possible eigenvalues are given by e„ — ±1. The latter 
fact is a consequence of the Eq. 

Pn = 1 (7) 

Letting e„ denote the parity of the n^^ plaquette, the 
9 possible states associated with the minor diagonal of a 
plaquette can be classified according to their parity, i.e. 
e„ = 1 for even parity states and e„ = — 1 for odd-parity 
states (see Table 1). 





State 


e 


2 holes 


|0> 


1 


bonding 


(cI..+4.J|0> 


1 


singlet 


Al,2|0) 


1 


antibonding 
triplet 


(cU-cL)lo) 

{clA^,, + s ^ s')\0) 


-1 
-1 



Table 1. Classification of the states of the minor 
diagonal of a plaquette according to their parity. 

A| 2 = (^c\ I + C2 1^ /\/2 is the pair field 
operator. 



The Hilbert space Tincckiacc of the t — J model can be 
split into a direct sum of subspaces He classified by the 
parity of their plaquettes, e = {en}n=i' namely 



necklace 



= ®e He 



(8) 



Every subspace Ti^ is left invariant under the action of the 
t—J Hamiltonian (Q) which can therefore be projected 
into an "effective" Hamiltonian Htj{e). In the previous 
section we have already seen an example of this type 
of decoupling phenomena. Indeed the alternating spin- 
1 /spin-i chain corresponds precisely to the case where all 
the plaquettes are odd and there are no holes. If holes 
are allowed then one has to consider, in addition to the 
triplets, the antibonding states on the odd plaquettes. 
Hence there are a total of 5 states at each site of the 
"effective" alternating chain associated with the minor 
diagonal of the odd plaquettes. 

On the other hand, if the parity of the plaquette is 
even, then the corresponding site in the chain has 4 pos- 
sible states which can be put into one-to-one correspon- 
dence with those of a Hubbard model as follows. 



Even diagonal 
empty(2 holes) 
bonding states (t, i) 
singlet state 



Hubbard site 
cmpty(l hole) 



single occupied state (t, J,) 
doubly occupied state 



(9) 



4 



In this case the Hamiltonian of the model is a t — J 
Hamiltonian with hopping parameter equal to \plt and 
with the same exchange and density-density couplings. 
There is no Hubbard U term. 

In summary, theorem (|^) implies that the necklace lad- 
der is in fact equivalent to a huge collection of alternating 
chains models where half of the sites are t — J like while 
the other half may be either spin 1 or spin-i antibonding 
states for odd parity, or Hubbard like, for even parity. 

It is beyond the scope of the present paper to study 
such an amazing variety of chain models disguised in the 
innocent looking necklace ladder. Instead, a more rea- 
sonable strategy is to ask for the values of the total spin 
S and parity e which give the absolute minimum of the 
GS energy, keeping fixed the values of the number of pla- 
quettes iV, the number of holes h and the ratio J/i, i.e. 



,J/t)<EQ{N,h,S,e,J/t), 



ys,e, 

(10) 



Even this question is not easy to answer with full gener- 
ality. However we shall study a few cases which suggests 
a general pattern for the behavior of the spin and parity 
as functions of doping. 



V. DMRG AND RVA RESULTS ON THE 
NECKLACE LADDER 

We shall concentrate on the case of a necklace ladder 
with 7 plaqucttcs and open BCs, which will allow us to 
present in a simple manner the basic features of the GS 
for various spins (S) and dopings (h). The values of the 
coupling constants are fixed to t = 1 and J — 0.35. 

In Table 2 we show, for several pairs {h, S) the parities 
of the plaquettes, the total GS energy computed with the 
DMRG and the RVA methods. This table also lists the 
label of the corresponding figures showing DMRG results 
for the hole and spin densities of the corresponding state. 
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FIG. 7. Pictorical representation of the most probable state 
for doped x = 1/3 necklace ladder of dimension 8x8. Blank 
circles denote holes and vertical solid lines represent valence 
bond states (Case (8,0)). 

The RVA results have been derived from an inhomoge- 
neous recurrence variational ansatz (see Appendix A). As 
in the latter cases we start from a state, hereafter called 
"classical" , which is considered to be the most important 
configuration present in the actual GS. Next we include 
the local quantum fluctuations around the classical state. 
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FIG. 8. Electronic densities for a necklace ladder as de- 
picted in Fig. 7. It corresponds to 7 plaquettes and 8 holes 
(doping X = 1/3). Below are shown the DMRG and RVA val- 
ues for the sites in the main diagonal, while above are shown 
the values on the rungs corresponding to subdiagonals. 

This is done for a whole set of "classical" states having 
the same number of plaquettes, holes and z-component 
of spin. As discussed in the appendix, the classification 
of the "classical" states is achieved by means of paths in a 
BratcUi diagram generated by folding and repeating the 
Dynkin diagram of the exceptional group Lie Eq. The 
six points of Eg are in one-to-one correspondence with 6 
different states on the necklace ladder, while the links of 
Eq are nearest neighbor constraints derived on the basis 
of the DMRG results in the region < a; < 1/3. The 
Bratelli construction gives us a systematic way to explore 
the GS manifold in the underdoped region x < 1/3. The 
overdoped region has to be studied with delocalized RVA 
states as discussed in references p^ , pjt . 
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(ei, • • • ,£7) 


Figure 
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(+ + + + + + +) 


7 and 8 


-16.554153 


-16.33996 


8 


1 


(-1- + H h + +) 


9 


-16.284855 


-16.00358 


7 


1/2 


(-f + + + + + -) 


10 


-15.489511 


-15.18141 


6 





(- + + + + + -) 


11 


-14.424805 


-14.02286 


6 


1 


(- + + + + + -) 


12 


-14.424798 


-14.02286 


9 


1/2 


(+ + + + + + +) 


13 


-16.746112 




10 





(+ + + + + + +) 


14 


-16.927899 




10 


1 


(+ + + + + + +) 


15 


-16.718476 





Table 2. DMRG and RVA total energies for a 
7-plaquette necklace ladder with h holes and total spin 
Sz- The string of epsilons is the pattern of parities in 
the subdiagonals (rungs). 
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Let us comment on the DMRG results shown in the 
figures. 

• Case (8,0): This is the most interesting case and 
it corresponds to one hole per site along the princi- 
pal diagonal. In the infinite length limit this state 
has doping x — 1/3. For this reason we shall here- 
after call this state the x — 1/3 state. Figure 7 
shows the most probable configuration which oc- 
curs when the holes occupy the principal diagonal 
of the ladder and the spins form perfect singlets 
along the minor diagonals. The latter fact implies 
that all the plaquettes are even (see Table 1). Fig- 
ure 8 shows the electronic density along the ladder 
computed with the DMRG and the RVA. 




Cases (6,0) and (6,1): (See Figs. 11-12) The 
state X = 1/3 is now doped with two electrons, 
which go to the boundary plaquettes which change 
their parity. There seems to be a small effective 
coupling between the two spin 1 /2 at the ends of the 
ladder, which lead to a breaking of the degeneracy 
between the triplet and the singlet. This is remi- 
niscent of the effective spin 1/2 at the ends of the 
Haldane and AKLT open spin chains ||lj]. There 
also exists a weak effective coupling that breaks the 
four fold degeneracy of the open chains. 




FIG. 11. Results from DMRG showing the necklace state 
with 6 holes and spin Sz ~ 0. 



FIG. 9. Results using DMRG showing the necklace state 
with 8 holes and spin Sz = 1- The diameter of the circles are 
proportional to the hole density, and the length of the arrows 
arc porportional to (Sz), accoring to the scale in the box. 



Case (8,1): (See Fig. 9) The spin excitation of the 
X — 1/3 state is given by a spin 1 magnon strongly 
localized on an odd parity plaquette located at the 
center of the ladder. The other plaquettes remain 
even and spinless. The value of the spin gap is given 
by As = 0.27 (DMRG) and 0.32 (RVA). 




FIG. 12. Results from DMRG showing the necklace state 
with 6 holes and spin Sz ~ 1- 




FIG. 10. DMRG results for the hole and spin densities of 
the necklace state with 7 holes and spin Sz ~ 1/2. 

• Case (7,1/2): (See Fig. 10) This case is obtained 
by doping the a; = 1/3 state with an electron. The 
additional electron goes into either of the boundary 
plaquettes. The corresponding plaquette changes 
its parity to —1. 



• Case (9,1/2): (See Fig. 13) The x = 1/3 state 
is doped with one hole. The parity of the plaque- 
ttes remain unchanged and the extra spin 1/2 delo- 
calizes along the whole system perhaps with some 
SDW component. The differences between this hole 
doped case and the electron doped case (7,1/2) are 
quite striking. 

• Case (10,0) (See Fig.l4) This state looks very 
much the same as the a; = 1/3 but just with more 
holes. 

• Case (10,1) (See Fig. 15) Same pattern as in the 
(9,1/2) case with the spin delocalized over the 
whole system. 

In summary the DMRG results clearly suggest the ex- 
istence of two distinct regimes corresponding to dopings 
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< X < 1/3 and x > 1/3. In the overdoped regime 
the plaquettes are always even while in the underdoped 
regime they can be even or odd. Phase separation into 
even and odd plaquettes may also be possible. Our re- 
sults at this moment are ambiguous and further numer- 
ical work is required. The most important result is the 
peculiar structure of the a; = 1/3 state, which we shall 
study further in the next two sections. 




^ 0.05 

FIG. 13. Results from DMRG showing the necklace state 
with 9 holes and spin Sz = 1/2. 



and Hartree-Fock Q calculations in large lattices and 
experimentally in some nickelates compounds. 

On the other hand the x — 1/3 state is a kind of ID 
generalization of the OS of 2 holes and 2 electrons on 
the 2x2 cluster discussed in reference [|l^ , in connection 
with the binding of holes in the 2-leg and higher-leg lad- 
ders. One can also use this local structure to build up a 
variational state of the 2-leg ladder valid for any doping 

The GS of two holes in a plaquette is the localized 
Cooper pair depicted in Fig. 16 and can be generated by 
the pair field operators acting on the vacuum as, 

I Cooper Pair >= (11) 
(a(A1_3 -t- Al, + + Al,) + {Al, + Al,)) |0) 

where A is given by 



[2+(J/4i)2]i/2- J/4t 



^ 0.0001 
O0.5 



FIG. 14. Results from DMRG showing the necklace state 
with 10 holes and spin Sz ~ 0. 




FIG. 15. Results from DMRG showing the necklace state 
with 10 holes and spin Sz — 1. 



VI. THE a; = 1/3 STATE OF THE NECKLACE 
LADDER 

The most important configuration contained in the x = 
1/3 state has spin singlets along the diagonal. This is 
consistent with and helps explain the tt phase shift in 
the (1,1) domain walls observed numerically with DMRG 



ICooper Pair > = a o°y + + s^o + o 

+ b fo I o + — ^ ) 
^1 o ^ 

FIG. 16. A pictorical representation of a localized Cooper 
pair on a plaquette. The parameter A in (|l^) is 
A = a/V2,b = 1. 

If J/t < 2, then A < 1 which means that a diagonal 
bond is more probable that a non-diagonal one. This 
feature is also observed in the a; = 1/3 state where the 
most probable bonds are those that line up along the 
transverse diagonals of the ladder ( Fig. 7) . Taking into 
account this state together with some important local 
fluctuations around it lead us to propose an RVA state 
with X = 1/3, which is defined by the recursion relations 
depicted in Fig. 17 (see Appendix A for details). The 
RVA state corresponding to = 1 plaquette coincides 
precisely with state (|l]) identifying A = a/ v^, 6=1 (See 
Fig. 16). For ladders with more than one plaquette {N > 
1) the symmetry of the diagonals disappear and 6 ^ 1 in 
general. In this case we have two independent variational 
parameters a and b. In Fig. 18 we give the energy per 
plaquette and the values of a and b as functions of the 
number of plaquettes of the necklace ladder obtained by 
minimization of the energy of the RVA state. We observe 
that both a and b become less than one in agreement with 
the DMRG results. All these results are quite satisfactory 
but still they do not give us a transparent physical picture 
of the x = 1/3 state. This will be done in the next 
section. 
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FIG. 17. A pictorical representation of the Recurrence Re- 
lations employed with the RVA method (See Appendix A) 
to construct variational G.S. states for the doped x = 1/3 
necklace ladder. The diagonal squares represent bulk states 
of a given length. Small circles represent holes and solid lines 
valence bonds. 



2-leg spin ladder can be seen in the strong coupling limit 
as the energy cost for breaking a bond along the rungs. 

In Section II we suggested that diagonal ladders could 
be thought of as collections of coupled plaquettes (Fig. 
3). The trouble is that in doing so one actually needs 
more sites than those available in the original lattice. 
Indeed the necklace ladder with N plaquettes has 3A^ 
sites for N large while the extended or decorated ladder 
shown on the right hand side of Fig. 3 has AN . 

The solution of this problem is achieved on physical 
grounds by defining on the lattice an extended t — J 
Hamiltonian which, in a certain strong coupling limit, 
becomes equivalent to the standard t — J Hamiltonian 
on . The extended Hamiltonian can also be studied 
in the limit where the plaquettes are weakly coupled. As 
we shall see, the latter limit provides a useful physical 
picture of the properties of the necklace ladder for x = 
1/3 and other dopings as well. 



A. The tp — Jp — td model 
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FIG. 18. Ground state energy per site as function of the 
number of plaquettes for a doped x = 1/2, necklace ladder. It 
is obtained using the RVA method (See Appendix A) . Here a 
and b denote the variational parameters and it is shown how 
they get stabilized towards their thermodynamic values as the 
length of the ladder increases. 



VII. THE PLAQUETTE PICTURE OF THE 
NECKLACE LADDER: THE tp - Jp - ta MODEL 

An interesting property of the rectangular ladders is 
that the strong coupling picture of the GS and excited 
states is generally valid also in the intermediate and weak 
coupling regimes. Thus, for example, the spin gap of the 



We shall define on the lattice an extended t 
model by the following Hamiltonian, 



J 



Hpd — 



Hp + Hd 



(13) 



= ^n,n+l(^d) 

where hn(tp, Jp) is a standard tp — Jp Hamiltonian in- 
volving only the 4 sites of the plaquette labelled by n. 
Of course /i„ and commute for n ^ m. On the 
other hand, hn,n+iitd) is a hopping Hamiltonian associ- 
ated with the link that connects the two nearest neighbor 
plaquettes n and n -I- 1 . Denoting by L and R the corre- 
sponding sites on the different plaquettes joined by the 
link < L, R > then hn,n+i is given by the link Hamilto- 
nian defined as 



td{c 



t 

L,s 



(14) 



B. Strong hopping limit of the tp — Jp — td model 

We want to prove that in the strong hopping limit, 
where — > oo, the tp~ Jp — td model becomes equivalent 
to the t — J model on the necklace ladder £^ . 

In this limit we first diagonalize Hd looking for the 
low energy modes of the plaquettes. We then define a 
renormalization group (RG) operator T, that leads to 
a renormalization of operators in the extended lattice 
model into operators that act on the necklace lattice. In 
particular the Hpd Hamiltonian is truncated to an ef- 
fective Hamiltonian which is equivalent to the original 
necklace Hamiltonian. The truncation operation is given 



8 



by the Eq. (for a review of the Real Space RG method 
see pOf) 



ff 



(15) 



The Hamiltonian hn.n+i{td) acts in a Hilbert space of 
dimension 3x3 = 9. It has two eigenvalues E = 
and 2td, with degeneracies 3 and 6 respectively. The 
zero eigenvalue corresponds to the states with two holes 
and the bonding state with up and down spins. In the 
limit td >> tp,Jp one retains only the latter 3 degrees 
of freedom which can be thought of as renormalized hole 
and spin up and spin down electron states, respectively. 
The truncation operator T, that maps the Hilbert space 
Hcp into the effective Hilbert space Hc^ is, given by 



T : Ti^p 

l°,-> 

|o,o) 





Til*) 
\o) 



(16) 



where • and o stand for one electron, with spin up or 
down, and one hole respectively living on a given link of 

, while * and o are the effective electron and holes liv- 
ing on the corresponding site of obtained by contract- 
ing the previous link to a site. The hermitean operator 

acts as follows, 



Tt : n, 



lo.o) 



(17) 



Eq.(p^ means that an electron state |*) of becomes 
the bonding state in the enlarged Hilbert space Ti-cp- 
The RG operators T and T'^ defined above satisfy the 
following Eqs. pi 



TT^ = 1 



(18) 



where Pq is a Gutzwiller operator which now acts on 
links rather than on sites as follows, 



while for i = M they act at the effective "middle" point 
of the link (i.e. < L,R >^ M). Of course, the opera- 
tors and states that are not on the principal diagonal of 
both and are not affected by the renormalization 
procedure. 

Using Eqs (|l5|) and (|^) we can immediately find that 
the renormalized effective Hamiltonian is given by the 
t — J Hamiltonian (^, i.e. 



-ffcff — T Hpd — T Hp{tp, Jp) — Hi 



(21) 



with the following values for the coupling constants. 



1 



V2 



J — ^ Jp, (necklace ladder) (22) 



In the derivation of (|2^) we are assuming periodic bound- 
ary conditions, along the principal diagonal of the neck- 
lace ladder. 

The strong hopping limit studied above is reminis- 
cent of the strong coupling limit of the Hubbard model 
which leads to the t-J model plus some extra three-site 
terms which are usually ignored. In the latter case the 
strong Coulomb repulsion forces the Gutzwiller on-site 
constraint. Our case is a "dual version" of this mech- 
anism, in the sense that the coupling constant involved 
is a hopping parameter, and that the Gutzwiller con- 
straint arises from a link rather than from a site con- 
straint. In the case of the tp — Jp — td model one does 
not have to do perturbation theory in order to produce 
the exchange term in the effective Hamiltonian since it 
is already contained in the plaquette Hamiltonian. Per- 
turbation theory would produce terms of the order 1/td, 
but they vanish at td — 00. The construction we have 
performed in this section can in principle be generalized 
to the Hubbard model |21|. 



The analogy between the Hubbard and the tp — Jp — td 
model suggests that we may learn something about the 
strong hopping limit by studying the weak hopping one. 
This is certainly true if there are no phase transitions 
between the two regimes. 



l°,-> 

lo.o) 



TLcp 


71 (I' 
71 (I' 

lo.o) 



(19) 



Using the above definitions we can easily obtain the 
renormalization of the different operators acting in , 



T CL,s T^^T CR^, Tt 
TulT^ =TnR Tt 



- -mCM,, 



1 



(20) 



Here Ci^s, Si, Ui are the fermion, spin and number opera- 
tors acting at the edges of the link < L,R> for i = L,R 



C. Weak Hopping Limit of the tp — Jp — td model 

In the weak hopping regime, i.e. td << tp, Jp, we first 
diagonalize the plaquette Hamiltonian Hp and treat Hd 
as a perturbation. The energy levels of Hp are given, to 
lowest order in perturbation theory, as tensor products 
of the eigenstates of every plaquette. There will be in 
general a huge degeneracy, which will be broken by the 
effective Hamiltonian derived from Hd using perturba- 
tion theory. Before going further into the study of the 
plaquette Hamiltonian we have to consider the relation- 
ship between the filling factors of the states belonging to 
lattices with different number of sites. 
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Let us consider a state in L with Nh holes and N^. 
electrons. Applying the operator T\ this state is trans- 
formed to a state in with 7V^^^ holes and ivj^^ elec- 
trons given by, 

nI^'^ = \nu + \n,, N,^p^=N, (23) 

These equations reflect the fact one gets an extra hole 
upon going to the enlarged lattice. Eqs.(p3|) imply 
the following relations between the doping factors x — 

Nh/iNh + iVe) and Xp = N^^p^ /{N^^p^ -f N.^p^), 

a; = i(4a;p-l), Xp = i(l + 3x) (24) 

From (^J) we get the following correspondences 

x = < — > Xp = 1/4 (25) 
X = 1/3 < — > Xp = 1/2 

which we shall discuss in detail below. 

Weak hopping picture of the x — 1/3 state 
Eq. (HI) implies that the 1/3 doped state of the neck- 
lace ladder is transformed to a state with two holes and 
two electrons per plaquette in the expanded necklace lat- 
tice. We show in appendix B that the lowest GS for this 
filling is given by the coherent superposition of Cooper 
pairs localized on the plaquettes, i.e. 

N 

\xp = 1/2, id = >= ]^ [Cooper pair >„ (26) 

where | Cooper pair >„ is the state given in Eq. . (pd]) with 
the parameters a and b given by Eq.([l^) for the values 
of tp^ Jp- 



the bonding and empty links will remain low in energy. 
In the limit when td becomes infinite we expect the GS 
( p6| ) to evolve continuously into the a: = 1/3 GS of the 
necklace. This suggest that the a: = 1/3 state of the neck- 
lace ladder can be described as a Gutzwiller projected 
state, i.e. 

\x = 1/3 T V^a'' = 1/2, td^O> (27) 

where we first project out the doubly occupied and an- 
tibonding states on the links on the principal diagonal 
of the expanded ladder and then project the resulting 
state into the Hilbert space of the necklace ladder. Some 
diagrammatics (Fig. 19) shows that the state (|2^ ) is ba- 
sically the same as the x = 1/3 RVA state constructed 
in Section VL This leads us to the conclusion that the 
X = 1/3 GS of the necklace ladder can be seen as the 
Gutzwiller projection of Cooper pairs localized on the 
plaquettes. In this case the Cooper pairs are locked in a 
Mott insulating phase and there is an exponential decay 
of the pair field. 

Weak hopping picture of the x — state 
The GS of a plaquette with 1 hole and 3 electrons 
for Jp/tp = 0.5 has spin 1/2 and it belongs to the 2- 
dimensional irrep labelled by E of the symmetry group 
D4 (see Appendix B). These two states differ in their 
parity along the minor diagonal which can be even or 
odd. Both states can be thought of as bound states of 
a Cooper pair and one electron ( see Fig. 20). The four 
fold degeneracy on every plaquette is broken by td- The 
odd parity plaquettes are lower in energy than the even 
ones and the effective model is given by a ferromagnetic 
spin 1/2 chain, 

Hcs = -'oft ^ Sjf • S^^^i (28) 

n 
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FIG. 19. Ground state for a doped x — 1/3 necklace ladder 
(down) obtained as the projection of a Xp = 1/2 doped state 
in the decorated (dual) diagonal ladder (above). 



Turning td on, the state (|2^) will be perturbed mainly 
along the principal diagonal. The doubly occupied and 
antibonding links will become high energy states while 



ICooper Pair + e' > = a 



+ 



+ b 



(— — ) 
^ o • ^ 



FIG. 20. A pictorical representation of a bound state 
formed by a Cooper pair and one electron. Small blank circles 
represent holes, black circles represent electrons and solid lines 
are valence bonds. Here d, b are relative amplitudes taken as 
variational parameters. 

Here S$f^ is the overall spin 1/2 operator of the odd pla- 
quette and Jeff ~ —t^ is a ferromagnetic exchange cou- 
pling constant. 

Following a reasoning similar to that for the case x = 
1/3 we conjecture that the x = state can be represented 
as the following Gutzwiller projected state, 

Ix = >~ T r^Q^ \xp = 1/4, e„ = -1, id = > (29) 



10 



whose structure is indeed very similar to the ferrimag- 
netic RVA state proposed in Section III. See Fig. 21 for 
a plaquette construction of the Neel state of the necklace 
ladder starting from the Xp — e„ = — 1 state. The 

gapless excitations of the ferrimagnetic a; = GS cor- 
respond, in the weak coupling picture, to the magnons 
of the ferromagnetic chain (^), while the gapped exci- 
tations correspond to an excitation of the plaquette to a 
state with spin 3/2. 



o 



0|| 



l(l) 



FIG. 21. A Neel- like ground state (down) for an undoped 
a; = necklace ladder obtained as the projection (above) of a 
doped state in the decorated (dual) diagonal ladder. 



A. The plaquette construction and medial graphs 

The plaquette construction of the necklace ladder is 
related to the so called medial graphs used in the color- 
ing problem or in Statistical Mechanics Before we 
show this connection we need to generalize our plaque- 
tte construction to diagonal ladders with more than one 
plaquette per unit cell , i.e. Up > 1. 

In this section we shall use the following notations. 



'-Up 



rectangular ladder with 
diagonal ladder with n 
(4, 8) lattice with Up plaquettes 



ng legs 
plaquettes 



(30) 



As an example we depict in Fig. 22 the lattices C2,C^ 
and /If . The lattice consist of 4-gons, i.e. plaquettes, 
joined by links, which are associated with the hopping 
parameter td while the plaquettes are associated with 
the parameters tp and Jp. For Up > 1 contains also 
8-gons that are formed by four td links and four tp links. 




In summary we have been able to obtain a satisfactory 
picture of both the a: = and 1/3 states in the weak 
coupling limit of the extended t ~ J model, which leads 
us to conclude that for these dopings there are no phase 
transition between the weak and strong coupling regimes. 
Other dopings involve the competition of the two elemen- 
tary plaquettes states used above and will be considered 
elsewhere. 




VIII. FROM ID TO 2D THROUGH DIAGONAL 
LADDERS 

The necklace ladder represents the first step in the di- 
agonal route to the 2D square lattice. In this section we 
shall push forward this viewpoint trying to see how much 
one can expect from it. This will lead us to ask questions 
whose solution we do not yet know. In this sense some 
of the material presented below is conjectural. 

Let us first start with a short excursion into graph 
theory. 



R 

L, = 



FIG. 22. Several examples of related ladders as explained 
in text. (Above) An example of (4, 8)-lattice with Wp = 2 pla- 
quettes. (Middle) An example of diagonal ladder with Up — 3 
plaquettes. (Down) An example of rectangular ladder with 
n=2 legs. 

As shown in the previous section the limit td ^ oo has 
the geometric significance of shrinking the corresponding 
trf-links into sites, so that the the lattice >C^^ "renormal- 
izes" into the diagonal ladder _i (see Fig. 22). In 
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this strong coupling limit the number of plaquettes actu- 
ally increases and some plaquettes are generated for free. 
The number of plaquettes of the diagonal ladder so ob- 
tained is odd. This construction does not produce even 
plaquette diagonal ladders. 

Observe that all the diagonal ladders are bipartite lat- 
tices but only when Up is even are the number of sites 
of the two different sublattices the same. This suggests 
that the rip-even diagonal AFH ladders belong to the 
same universality class as the n;-even rectangular lad- 
ders, while the n^-odd diagonal ladders belong to a dif- 
ferent universality class characterized by ferrimagnetic 
GS's. 

The opposite limit, where td ^ 0, has the geometrical 
meaning of shrinking the plaquettes into sites, so that 
"renormalizes" the system to a rectangular ladder 
with Up legs ( see Fig. 22). 

We summarize the above geometric RG operations in 
the following symbolic manner, 

< — At {u - 0) (31) 

In this sense the (4, 8) lattice is an interpolating structure 
between diagonal and rectangular lattices. 

There is an interesting connection between this plaque- 
tte construction and the theory of medial graphs. Con- 
sider a graph G made of a set of points i connected by 
links < i,j >. A medial graph M{G), associated with 
the graph G, is obtained by surrounding every site i of 
G by a polygon Pi, such that two polygons Pi and Pj, 
which correspond to a link < i,j >, meet at a single in- 
tersection point Pi n Pj , which lies on the middle of the 
link < i,j > ||2^ (see Fig. 23 for a generic example). 

Choosing the polygons Pi to be 4-gons, i.e. plaquettes, 
one can easily show that a diagonal ladder with an odd 
number of plaquettes is the medial graph of a rectangular 
ladder, namely 

£f„^_i= MiCl) (32) 

Medial graphs are used in Statistical Mechanics to 
show the equivalence between the Potts model and the 
6- vertex model |2jj2^. Indeed one can show that the 
Potts model defined on a graph G is equivalent, i.e. has 
the same partition function after appropriate identifica- 
tion of parameters, to the 6 vertex model defined on the 
medial graph A4{G), i.e. 

^Potts(G) = Zq 

—vertex 

{M{G)) (33) 

The transformation G M{G) is a kind of duality map 
that relates two seemingly unrelated models and it is in 
fact the key to solve the 2D critical Potts model in terms 
of the 6 vertex one. 




FIG. 23. An example of medial graph construction. Here 
the graph G is made by solid lines and blank circles. Its 
associated medial graph A4{G) is made by dashed lines. 



B. Plaquette Construction of the 2D Square Lattice 

In Fig. 24 we apply the plaquette construction to 
the 2D square lattice. It is a simple generalization of 
the construction shown in the previous subsection when 
rip — > oo. If £^ is a square lattice with lattice spacing 
a, then is also a square lattice but with spacing \/2a 



FIG. 24. Plaquette constrution (small interior squares plus 
dashed lines) of the 2D square lattice as explained in text. 

Let X be the doping of a t — J model defined on , 
and Xp the doping factor of a, tp — Jp — ta model defined 
on then the relations between these quantities are 
analogous to Eqs.(p^) and ( p^ for the necklace ladder, 
namely 

X = (2xp - 1) (34) 
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Eq.(Q) implies that the undoped system x = corre- 
sponds to doping Xp = 1/2 in the enlarged lattice. Fig. 
25 shows a plaquette construction of the Neel state from 
the Xp = 1/2 state. Notice that the plaquettes have spin 
1 and that the parity on their diagonals alternate be- 
tween (1,-1) and (—1,1). In the strong hopping limit 
the plaquettes have an effective spin 1. The whole set of 
these effective spin I's are coupled antiferromagnetically 
and form a square lattice with lattice spacing which is \/2 
times larger than the lattice spacing of the original spin 
1/2 model. In a certain sense the plaquette construction 
integrates out degrees of freedom and renormalizes the 
system into an AF Heisenberg model with spin 1 and 
lattice space ^/2a. This picture agrees qualitatively with 
the RG flow of the 0(3) non linear sigma model in the 
renormalized classical region at zero temperature |20] . 
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FIG. 25. Plaquette constrution of the 2D Neel state in 
a square lattice as explained in text. The corresponding 
depicted state on the decorated (dual) lattice has doping 
Xp = 1/2 and spin 1 on every small square plaquette. 

In the weak hopping limit, however the GS of the 2D 
model is given essentially by the coherent superposition 
of localized Cooper pairs used in the construction of the 
X = 1/3 necklace state. The Gutzwiller projection of this 
state onto the original lattice will produce Spin Peierls 
state rather than an AFLRO state. 

We conclude that unlike the case of the necklace ladder, 
the tp — Jp — td model in 2D must have a phase transi- 
tion for some intermediate value of td- The study of this 
model may serve to clarify the relationship between the 
AFLRO and the d-wave pairing structures observed in 
the theoretical models of strongly correlated systems. 



IX. CONCLUSIONS 

Diagonal t — J ladders provide an alternative route of 
interpolating between 1 and 2 spatial dimensions. Here 
we have described a general framework for such an inter- 
polation and introduced a generalized tp — Jp — td plaque- 
tte model in which the individual plaquettes are linked 
by a hopping term td- In the strong hopping td ^ tp, Jp 
limit, the generalized plaquette model was shown to map 
into the original diagonal t — J model with renormal- 
ized parameters and filling factor. Thus, the generalized 
tp — Jp — td model provides a dual model to the original 
diagonal model. In this sense, it is interesting to study 
the tp — Jp — td model in the weak hopping limit. If 
there is no phase transition between the weak and strong 
hopping limits, then the weak hopping limit can provide 
new insight into the nature of the original t — J diagonal 
ladder. We believe that this is the case for the Up — 1 
plaquette necklace ladder and that its groundstate for a 
doping X = ^ can be understood as the Gutzwiller pro- 
jection of a product state of Cooper pairs localized on 
the plaquettes of the quarter-filled extended tp — Jp — td 
model. Alternatively, for the rip ^ oo 2D limit, we be- 
lieve that the extended tp — Jp — td model at a doping 
Xp = ^, which corresponds to the undoped {x = 0) t — J 
model, will have a phase transition for an intermediate 
value of td/tp. In this case, our conjecture is that the 
strong coupling limit will have a ground state with long 
range AF order while the weak coupling phase will be a 
localized Spin Peierls state. 

In order to make these ideas more concrete, we have 
focused on the single plaquette rip = 1 necklace ladder. 
Here, using the results of DMRG and RVA calculations, 
we have studied the necklace ladder for various dopings x. 
For X = ^, the DMRG calculations show that in the most 
probable configuration, the holes occupy the sites along 
the principal diagonal of the necklace and the spins form 
perfect singlets along the minor diagonals. The RVA cal- 
culations, starting from a "classical" configuration and 
mixing in local quantum fluctuations about this state, 
provide a ground state energy in good agreement with 
the DMRG result. Then, as discussed above, a more 
transparent physical picture of the x = ^ state of the di- 
agonal necklace is provided by the extended tp — Jp — td 
dual model at a filling of which this state is 

seen as a localized Cooper pair state. It will be interest- 
ing to understand what happens when additional holes 
are added. In particular, will a necklace with a doping of 
Xp = 0.5 + S have power law d-wave like pairing correla- 
tions? 

For X — 0, the diagonal necklace is equivalent to an al- 
ternating s = 1/s = i spin chain and has a ferrimagnetic 
ground state with total spin N/2, with N the number 
of unit cells of the necklace. There are gapped excita- 
tions with spin N/2 + 1 and gapless excitations with spin 
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N/2 — 1. In the weak hopping hmit of the tp — Jp — td 
model, these excitations correspond to local excitations 
of the plaquettes to a spin 3/2 state and to magnons of 
a ferromagnetic spin ^ chain respectively. Thus, in the 
a; = case, the dual model provides a useful physical 
picture. 

We also have found that when the x — ^ state is doped 
with holes, the ground state plaquettes retain the even 
parity characteristic of the x = ^ state. However, when 
electrons are added, this parity can be even or odd. Thus, 
it appears that the x = ^ doping separates the system 
into two distinct regions. 

Clearly, the diagonal ladders form a rich class of mod- 
els with properties ranging from ferrimagnetic to anti- 
ferromagnetic and from localized pair states to possible 
extended pairing states. Furthermore, the tp — Jp — td 
model provides a dual description which suggests alterna- 
tive physical pictures and approximation schemes as well 
as connections to concepts from statistical mechanics. 
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APPENDIX A: RVA APPROACH TO THE 
NECKLACE LADDER 

The RVA method is a kind of simplified DMRG where 
one retains a single state as the best candidate for the 
GS of the system. As in the DMRG the GS of a given 
length is constructed recursively from the states defined 
in previous steps. This idea can be implemented analyt- 
ically if the ansatz is sufficiently simple. Below we shall 
propose various RVA states for the necklace ladder with 
dopings < a; < 1/3. 



1 3 5 2N-1 2N+1 




2N 



FIG. 26. A pictorical view of the necklace ladder showing 
the labelling convention employed to denote the variational 
RVA states. Here the positions along the main (horizontal) 
diagonal of the ladder are sing odd sites, while the even posi- 
tions are made up of rungs or sudiagonals. 



A. Case x = 

Let us begin by labelling the sites of the necklace ladder 
as in Fig. 26. The even sites denote the minor diagonal of 
the ladder while the odd sites are those on the principal 
diagonal. 

At zero doping there are only two possible states on 
the odd sites given by, 

I T>= c] |0 >, I i>= 4 |0 > (36) 

On the even sites there are a triplet and a singlet state 
given by, 

\S>^ Al_2|0> 
in >= cl cl, 10 > 

The Neel state on the necklace ladder can be written 
simply as, 

U> \T^> \i> \T^> • • • I i> |Tt > I i> |Tt > 

(38) 

A trivial observation is that the Neel state on the necklace 
with N sites is generated by the first order recurrence 
relation, 

|Neel,2A + 1 I i> |Neel,2iV> (39) 
|Neel, 2N |T| > |Neel, 2A - 1 > 

Quantum fluctuations around the Neel state amount to 
local changes of the form (see Fig. 27) 

U> ir^ >— . |a,TT) >= I T> |To> 

I^T > I i>— > im,i) >= iTo > I T> (40) 

I i> > I liLT^A) >= I T> in > I T> 

Hence a globally perturbed Neel state is a coherent su- 
perposition of states of the form, 

i (n i)m i) tt (i i) (n m i m d i 

(41) 

where the parenthesis denote the quantum fluctuations 
given in (^) . The RVA state is a linear superposition of 
states of the type ( ^ ) weighted with probability ampli- 
tudes which are the variational parameters. 

The RVA state is generated by the following recursion 
relations (RR), 

\2N+1 >= I i> \2N > 

+u |(i, T|) > |27V - 1 > 
+v\{i,T^,i)>\2N~2> (42) 
|2iV >= m > |2A - 1 > 

+u|(TT,i) > ;|2A-2> 
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with the initial condition 

|7V-1>=U>, |iV: 



>= 1 



(43) 



To compute the energy of the state > we define the 
foUowing matrix elements, 



Zn N\N > 
En =< N\Hn\N > 



(44) 



where is the Hamiltonian of the system with N sites. 

The RR's for the states ( ^2| ) imply a set of recursion 
relations for the matrix elements ([44|). Form the norm 
we get, 



Z2N+I — Z2N + U^Z2N-l + v'^Z2N-2 



Z2N — Z2N-I + 

The initial conditions are, 
Zq — Zi 

The RR's for the energy are. 



(45) 



(46) 



E2N+I = E2N + V?E2N-1 + v'^E2N-2 + 

J{^/2u - \)Z2N-1 + J{2V2UV - V^)Z2N-2 
+ WZ2N-?. 



E2N — E2N-1 + u E2N-2+ 



Wz 



2N-\ 



J{\f2u - ^)Z2N-2 + Ju'^Z2N-3Y' 

while in this case the initial data are 

Eo = Ei = (47) 

Minimizing the GS energy in the limit iV > > 1 we find 
that the GS per plaquette is given by -0.4822 J, which 
corresponds to an energy per site of the associated spin 
chain equal to -0.7233 J. The values of the variational 
parameters are given by u = —0.3288 and v = 0.1691. 



B. Case x = l/3 

The most probable state for this doping is given by 
(see Fig. 7) 



> \S> |o> \S> •••|o> \S> |o> 



(48) 



Analogously to Eq.(Q) we define local fluctuations 
around ( ^ ) in terms of the states Koj-S*) >,|(S', o) > 
and |(o,S',o) > depicted in Fig. 27. The x = 1/3 RVA 
state can then be constructed from the following RR's ( 
see Fig. 17) 



|2iV+ 1 >= |o > \2N > 

+a |(o, S)> |2A^ - 1 > 
+b |(o,S',o) > |2iV - 2 > 
\2N >= \S> |2iV - 1 > 

+a 1(5, o) > |2iV-2 > 



(49) 



The norm of |iV > satisfies the RR's (|45]) with the 
replacements u a,v ^ b. The RR's for the energy 
matrix elements are given by, 

E2N+1 = E2N + a^E2N-i + b'^E2N-2 (50) 
+ (V2(-2t)a - Ja'^)Z2N-i + V2{-2t)2abZ2N-2 
+ a^{-J/4)(a^Z2N-3 + b^Z2N-4) 

+ &2(- J/4)(2Z2JV-3 + a^Z2N-4) 

E2N — E2N-1 + o?E2N-2 — (2\f2ta + Ja^)Z2N-2 

i-J/4)ia^Z2N-3 + 2b^Z2N-i) 
+ a^-J/4)(2Z2N-3 + a^Z2N-4) 

The initial conditions for both Zn and En are the same 
as for the undoped case. In the limit iV — > 00 we find 

lim ^2^^ = -0.7387, a = 0.7873, b = 0.5478 



(51) 



which give the asymptotic values of the curves in Fig. 18. 



I(0,S)> = +/ ^1/2 1(0, T+)> = (I I - f ^ ) 1/^ 



l(|,TV)> 



FIG. 27. A pictorical representation of fluctuation states 
as constructed in text for the RVA method in the necklace 
ladder. Here r+ and T- stand for Tj and T^, respectively. 



C. Cases < a; < 1/3 

In the underdoped region < x < 1/3 we have ob- 
served with the DMRG that many of the GS that one 
gets, and particularly those listed in Table 2 can be un- 
derstood as quantum fluctuations around a "classical" 
state IV'o >• This state has the generic structure already 
seen in the cases x — and 1/3 (see (^ ) and (^8|)), 
namely 



|V'0) = |/2JV+l>|M-- -1^2)1/1) 



(52) 



where the states i = 1, 2, • • • 2A^, 2A^ -I- 1, are taken 
to be 

|/odrf)={|0),U>, IT)} = {1,3,5} (53) 



|Uen)={|5),|TT),|rx)} = {2,4,6} (54) 
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Notice that we do not allow the holes on the minor diag- 
onals of the classical state \ipo >■ They will go there after 
considering the fluctuations [Q. Based on the DMRG 
results as well as physical considerations, we shall allow 
the following pairs \li > > in \tPq > 



ms), io)|tt), \m), \om), 



IT, 



|5)|0), |rT)|0), \T^)\l), \T,)\0), |T,)|T) 



(55) 



(56) 



l(S,0)> = ( / + \) 1/2 l(T+, 0)>= ( I + - ^ i ) lY5 
l(T.,0)> = (°^-i ^)., KT. ,!)>=( I 

i(T+4)>= (J ^ + 1 

FIG. 28. A pictorical representation of fluctuation states 
as constructed in text for the RVA method in the necklace 
ladder. Here T+ and T- stand for and Tj, respectively. 



O 



l(0,S,0)>: 



o 



l(0,T+,0)>= t f l(0,X,0)>=i i 

O o ' 

X i 

l(|,T+,|)>= t t l(t,T.,t)>=| I 

l(0,T4)>= I ^ ^ = l(|,T+,0)> 



1(0, T.,l)> = fl ° 1 - I \ )l/f2 = l(t 'T. , 
^1 O ^ 



0)> 



FIG. 29. A pictorical representation of fluctuation states 
as constructed in text for the RVA method in the necklace 
ladder. Here T+ and T_ stand for Tj and Ti, respectively. 

This connectivity of the states making up a certain \^q) 
state can be summarized in a graph in which we place a 
site for each and every 6 states in (53), (|5^), and joint 
them by links whenever it is possible to find them one 
next to each other in the lip^) state according to the al- 
lowed local configurations (^ 56). This graph is depicted 
in Fig. 30, and coincides with the Dynkin diagram of the 
exceptional Lie Group Eq. 



Now we can characterize every admissible classical 
state {ipo) in a geometrical fashion: each j-^o) is a path in 
the so called Bratelli diagram associated to the Dynkin 
diagram of Eq. This Bratelli diagram is shown in Fig. 
31. The way it is constructed is apparent in that figure: 
one starts with the 3 possible site-states (^) located one 
on top of each other. These states are located by the 
label li of the first site of the diagonal ladder. Then we 
link them to the 3 possible rung-states (|4|) according to 
the connectivity prescribed in Fig. 30. These rung-states 
are located by the label I2 of the second position of the 
diagonal ladder. Once this is achieved, the rest of the 
graph in Fig. 31 is made up by reflecting this basic piece 
over the rest of the labels ^3, Z4, . . . , /27V, hN+i- Observe 
that the x — and x = 1/3 states discussed previously 
correspond to straight paths of the Bratelli diagram (|2S|). 
A similar type of construction is also used in Statistical 
Mechanics in the context of the Face Models lESl . 



o- 
I \ > 



o — 



o 



-o- 



i> 

o — 



IO> 



-o 

It 



FIG. 30. Dynkin diagram of the exceptional Lie group Eg 
and its associated site and rung states contributing to the 
variational RVA method for the necklace ladder. These 6 
states make up every classical state \tpo) on the underdoped 
ladder according to the connectivity of this diagram. 

The quantum fluctuations around jV'o > amounts 
to considering the normalized states \{li,li+i) > and 
\{li,li+iJi+2) > depicted in Figs. 27-29. An interest- 
ing property of these states is that they are orthogonal. 



{lj\{li,k+i,li+2)) =0, j ^i,i + l,i + 2 



(57) 



{{lj,l,+i)\{k,k+i,h+2)) ^0, + l (58) 

The RVA state built upon I'lpo > is generated by the RR's, 
(see Fig. 30) 

\2N+1) - \l2N+i)\2N)+aN\il2N+ij2N))\2N-l) 

+bN\il2N+l,l2N, l2N-l))\2N - 2) 

|2iV) = |;2w)|2iV-I)+CAr|(;2iv,;2JV-i))|2iV-2) 
provided with the initial data. 



(59) 
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10)^1 



(60) 



7 (r? ?.2 JJ) 

^27V-4(,aArOjV-1^2Ar+l,2Ar,27V-1.2Af- 



2,27V-3 



Using the orthogonality conditions ( |57| , |58| ) it is easy to 
get the RR's satisfied by the norm of the RVA state 

^27V+l = Z2N + a^Z2Ar_i + bj^Z2N-2 (61) 



"N^N~l'^2N+l,2N.2N-l,2N-2.2N-3) 



E2N — E2N-I + Cj^E2N-2 



(63) 



^2N 



Z2N-I + C^Z2N-2 



(62) 



7 / (0) I 2 (2) . r, (1) \ 

-^2N-2((-2N,2N-l + '^N^2N,2N-1 + '^'^N^2N ,2N -1)' 



Z2N-3{aN-l('2N,2N-l,2N 




FIG. 31. Bratelli diagram of the exceptional Lie group Ed- 
it serves to classify all the classical states i/jq appearing in 
the RVA treatment of the underdoped necklace ladder: every 
path on this diagram characterizes one of those classical states 
and provides its quantum numbers. Here r+ and r_ stand 
for and T;, respectively. Also, S denotes a singlet state 
and O represents a hole. 

The RR's for the energy have the structure given below 

E2N+I = E2N + a%E2N-l + b%E2N-2 + 



^2N-\y^2N-\-l,2N ^ "W=2W+1,2W 



+2aive27V+i,2Ar)+ 



7 (r"^ f^"*) -X- 

^2W-2lC^e2JV+l,2JV,2Ar-l ' 



o2 f W , r.2 (8) , 

"A''*'2Af+l,2Af,2W-l ^ "N'^2N+\,2N,2N-\^ 



'^^NCNf'2N+l,2N,2N-l 2^^^^^2Af+l,2Af,2Af-l) 



7 / 2 2 (5) I 

^27V-3iaAr Oat- 1 £2^+1, 2Af,2Af-l,27V-2~r 



"N''2N+1,2N.2N-1,2N-2)^ 



~^'^N^2N,2N~l,2N-2)'^ 



Z2 



2N-i{b'N-1^2N,2N-l,2N-2.2N-3 



I 2 2 ^ 

+ ^N(^N-1^2N,2N-1,2N-2,2N-3) 



.(5) 



(64) 



The e symbols in Eqs. (|64| ) are reduced energy matrix 
elements involving the states defined in Fig. 27-29. To 
simplify the RR's for the energy we adopt the notation: 
e2N+i,2N,... = Q2iv+i,i2N,...i where I2N+1J2N, ■ ■ ■ mn over 
the six possible states, 1, . . . , 6, of the Dynkin diagram 
Eq . All the non vanishing values of e are shown in tables 
3. 






FIG. 32. A pictorical representation of the Recursion Rela- 
tions in Eq. (M) employed to generate the variational states 
in the RVA treatment of the underdoped (0 < a; < 1/3) neck- 
lace ladder. Here a, b and c are local variational parameters. 
A square denotes a bulk state on a ladder of a length given by 
its number inside. The black circles and solid lines represent 
generic fluctuation states as explained in the text. 

Given the previous equations we set up the following 
strategy to derive the results presented in Section V. 

• Fix the length N of the ladder, the number of holes 
h and the third component of the spin of the 
whole ladder. 
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• Generate all the \ipo > 

configurations with those quantum numbers 
N, h, . Generically the number of configurations 
grows exponentially. For example the 7 plaquette 
case studied in Section V has N=15 and a total of 
3 X 10^ configurations. 

• Compute the energy of the state associated with 
the zero-order state IV'o > using the recursion re- 
lations and find the variational parameters which 
lead to a minimum energy. For example for the 
7 plaquette ladder we used 21 independent varia- 
tional parameters. 

• Extract the state > which has the absolute min- 
imum energy for a given N, h and . 





(3,4), (5,6), (4,3), (6,5) 


-J 




(1,2),(1,4),(1,6),(2,1),(4,1),(6,1) 


-^/2^ 




(3,4),(5,6),(4,3),(6,5) 


J 

V2 


n .12 


(1,2),(2,1) 


-J 


n ,'i2 


(3,4),(5,6),(4,3),(6,5) 


J 
2 



Table 3 a) Non- vanishing reduced energy matrix 

elements for ef^,^^ = {li.2\{hi\h^i^f^\li^)\li^) where the tJ 
Hamiltonian acts only on the states specified by its 
subscripts. Likewise for the elements 



il.42,«3 


(1,2,1) 


-V2t 




(1,4,3),(1,6,5) 


-t 




(3,4,1),(5,6,1) 


-,/ 

2 


e^^^. - 


(3,4,3),(5,6,5) 


J 




(1,2,1) 




'^il.i2-i3 


(1,4,3), (1,6,5) 


J 

2 


«l,«2,'i3 


(3,4,1), (5,6,1) 


t 




(3,4,3), (5,6,5) 


J 

V2 



Table 3 b) Non-vanishing reduced energy matrix 

elements for e^l]^^^^^ = {{ln,li2Jh)\h\l2\ln)\{k2,ks)) 

where the t J Hamiltonian acts only on the states 
specified by its subscripts. Likewise for the elements 

^lfL.„ = ((/a-',2-'M)i/'!!-!i(/u-',2)>i/.3>- 



«l,42,'i3 


(2,1,2), (2,1,4),(2,1,6),(4,1,2),(6,1,2) 


,7 
2 


41,42,13 


(3,4,3),(5,6,5),(3,4,1),(5,6,1) 


J 
2 




(6,1,4), (4,1,6) 


J 


^41,42,43 


(2,1,2),(2,1,4),(2,1,6),(4,1,2),(6,1,2) 


,7 
2 


41,42,43 


(3,4,3),(5,6,5),(1,4,3),(1,6,5) 


,7 
2 


4l,42,».'i 


(4,1,6),(6,1,4) 


J 



Table 3 c) Non-vanishing reduced energy matrix 

elements for e^^fj.^ = ((/i^, ^43)l(^n I^4i)|(i42, ^43)) 

where the t J Hamiltonian acts only on the states 
specified by its subscripts. Likewise for the elements 

t]i2M = iliMln,li2)\h^l\{li.,li2))\hs)- 



41,42,43 


(3,4,3),(5,6,5) 


-2J 


41,42,43 


(1,4,3),(1,6,5),(3,4,1),(5,6,1) 


-J 



Table 3 d) Non-vanishing energy matrix elements 

e£i2,43 = {ih.Ji2Jis)\{h'^J^+ht2)\{k,,k,,k,)) where 
the t J Hamiltonian acts only on the states specified by 
its subscripts. . 



41,42,43,44 


(1,2,1,2),(1,2,1,4),(1,2,1,6),(1,4,1,2),(1,6,1,2) 


J 

4 


,(^). . . 
41,42,43,44 


(1,4,1,6),(1,6,1,4),(3,4,3,4),(5,6,5,6),(3,4,1,2) 


J 
2 


41,42,43,44 


(5,6,1,2),(3,4,1,6),(5,6,1,4),(3,4,1,4),(5,6,1,6) 


,7 
2 


^(5). . . 
41,42,43,44 


(2,1,2,1),(2,1,4,1),(2,1,6,1),(4,1,2,1),(6,1,2,1) 


J 

4 


^(5). . . 
41,42.43,44 


(2,1,4,3),(2,1,6,5),(4,3,4,3),(6,5,6,5),(4,1,4,3) 


J 
2 


41,42,43,14 


(6,1,6,5),(4,1,6,1),(6,1,4,1),(4,1,6,5),(6,1,4,3) 


J 
2 


Table 3 e) Non-vanishing energy matrix elements 

^41,42,43,44 ~ ((^43, ^44)l((^4l, ^42)1^8243 1 ('h , ^42))l('43, ^u)) 

where the t J Hamiltonian acts only on the states 
specified by its subscripts. 


41,42.43,44 


(2,1,2,1),(2,1,4,1),(2,1,6,1),(2,1,4,3),(2,1,6,5) 


J 

2 


41,42,43,44 


(4,1,2,1),(6,1,2,1),(4,1,6,5),(6,1,4,3) 


J 
2 


,(6), 

41,42,43,44 


(4,1,6,1),(6,1,4,1) 


-J 


Table 3 f) Non-vanishing energy matrix elements 

4i!42,43,44 = (('42,k3,^44)l(^4i|/i!:5|/4i)|(/42,k3,^44)) whCIC 

the tJ Hamiltonian acts only on the states specified by 
its subscripts. 


41,42,43,44 


(1,2,1,2),(1,2,1,4),(1,2,1,6),(1,4,1,2),(1,6,1,2) 


,7 
2 


41,42,43,44 


(3,4,1,2),(5,6,1,2),(3,4,1,6),(5,6,1,4) 


J 
2 


41,42,43,44 


(1,4,1,6), (1,6, 1,4) 


-J 



Table 3 g) Non- vanishing energy matrix elements 

^ii!42,43,44 = (^44l(('4l,^42,^43)|/l!3/il('4i,^42,^43))|'44) whcrC 

the t J Hamiltonian acts only on the states specified by 
its subscripts. 



e^'l . . . 

tl ,l2 -^■^ ,«4 


(1,2,1,2,1),(1,2,1,4,1),(1,2,1,6,1),(1,2,1,4,3) 


J 

4 


. . . 

41,42,43,44,45 


(1,2,1,6,5),(1,4,1,2,1),(1,6,1,2,1) 


,7 
4 


41,42,43,44,45 


(1,4,1,6,1),(1,6,1,4,1),(1,4,1,6,5),(1,6,1,4,3) 


,7 
2 


(7) 

41,42,43,44,45 


(3,4,3,4,3),(5,6,5,6,5),(3,4,3,4,1),(5,6,5,6,1) 


J 
2 


41,42,43,44,45 


(3,4,1,2,1), (5,6,1,2,1),(3,4,1,4,3), (5,6,1,6,5) 


.7 
2 


41,42,43,44,45 


(3,4,1,4,1), (5,6,1,6,1), (3,4,1,6,5) 


.7 
2 


,(7) 

41,42,43,44,45 


(5,6,1,4,3), (3,4,1,6,1), (5,6,1,4,1) 


J 
2 



Table 3 h) Non- vanishing reduced energy matrix 



elements for e 



(7) 



iiks , hi Jt,)\ {{hi , k2 I (^^1 ' )) I (ha , hi ,k,)) where 

the tJ Hamiltonian acts only on the states specified by 
its subscripts. 



11,22, 43. 14, «5 


(1,2,1,2,1), (1,4,1,2,1),(1,6,1,2,1), (3,4,1,2,1) 


.7 
4 


. . . 


(5,6,1,2,1), (1,2,1,4,1),(1,2,1,6,1) 


,7 
4 


. . . 

Jl,i2,i3,U,«5 


(1,6,1,4,1), (1,4,1,6,1),(5,6,1,4,1), (3,4,1,6,1) 


,7 
2 


e-(^) . . . 

Il,l2,l3,«4,l5 


(3,4,3,4,3), (5,6,5,6,5),(1,4,3,4,3), (1,6,5,6,5) 


,7 
2 


Jl,i2,«3,l4,l5 


(1,2,1,4,3),(1,2,1,6,5),(3,4,1,4,3), (5,6,1,6,5) 


J 
2 


e-(^) . . . 

Il,*2,*3,«4,»5 


(1,4,1,4,3), (1,6,1,6,5),(5,6,1,4,3) 


J 
2 


,-(7) 


(3,4,1,6,5), (1,6,1,4,3),(1,4,1,6,5) 


.7 
2 



Table 3 i) Non-vanishing reduced energy matrix 



elements for e,^ , - - - 



{{K iih^ > ^2 , ^23 ) I hf/il \ilh,k2,h3))\ {K ^hrj) where 
the tJ Hamiltonian acts only on the states specified by 
its subscripts. 



D. Electronic Density G.S. Expectation Values 

Here we derive the general equations to compute 
the electronic density expectation values using the RVA 
method. In particular, we have plotted in Fig. 8 the re- 
sults for the doping case x = 1/3 in the necklace ladder 
and compared it with the DMRG results showing a good 
agreement. 

Let us denote the density electronic operator at the 
site position R as dfl, namely. 



Using the RR's for the diagonal (|9|) ladder we may 
find also RR's for the unnormalized V.E.V.'s: 

d2N+liR) - d^NiR) + Ci%d2N-l{R) + b%d2N-2{R) (68) 



d2NiR) = d2N-liR) + C%d2N-2iR) 



(69) 



whenerver the position R of the insertion is not near the 
end of the diagonal ladder state. The derivation of these 
RR's follow closely that of the norms Z's using the or- 
thogonality relations. 

Now we need to determine the boundary or initial con- 
ditions to feed those RR's. Let us consider the cases R 
odd and even separately. 



1. R Odd 

As for the odd sites the RR's of the states are third 
order, we have to compute directly the values (i/j(i?), 
dR+i{R) and dji^2{R)- Clearly we also have, 



dR-i{R)=0 



(70) 



dR = Ur^^ + HR^l 



(65) 



Now, as R is odd, let us set R = 2r + 1. Then, using the 
RR's we find, 

d2r+l{2r + 1) = 4° + 1^2r + a'^2r-l + &'^2,.-2 (71) 

In the derivation of these RR's we need to compute the 
following reduced density matrix elements, 

_ /, I J I, \ - / if ^2r+l = |o> 

^2r+l = (^2r+l|d2.+ l|/2,.+ l) - | ^ if l^^^^^ = || | 

(72) 



'^2° + l,2r = ((^2r+l, ?2r)|d2r+l|(?2r+l, ?2r)) — 1 (73) 



where nR,ij and tir.i are the number operators for 
fermions with spins up and down, respectively. 
For a diagonal ladder of length N (sites-|-rungs) we need 
to compute the V.E.V. of the density operator dR in the 
ground state, which we denote as. 



dN{R) - {N\dR\N) 



(66) 



When the position R is even we shall need an additional 
index k = 1, 2 to locate the upper site (fc = 1) and the 
lower site {k = 2) of the rung. 

Notice that these density V.E.V. ( |66|) are not normalized. 
We may introduce normalized densities as. 



dN{R) 



{N\dR\N) ^ dN{R) 
{N\N) Zn 



(67) 



The densities dM{R) takes on values from to 1 de- 
pending on whether we find a hole with maximum prob- 
abihty (djv(^) — 0) or one electron {d^iR) — 1). 



"2r+l,2r,2r-l — 
{{hr+l,l2r,l2r-l)\d2r+l\{hr+l,hr,l2r-l)) = 1 (74) 

Likewise we proceed with the initial value dR+i (R) and 
we find, 

d2,.+2(2r + 1) ^ d2r+i(2r + 1) + c^+iZa. (75) 
where we have used the result, 

((^2r+2, ^2r+l)M2r-|-l|('2r+2, ^2r+l)) — 1 (76) 

Finally, we proceed similarly with the other initial 
value dR^2{R) and we find, 

d2r+3(2r + 1) = d2r+2(2r + 1) + 



al+id2r+i{2r + 1) + fo^+l^Zr 



(77) 
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2. R Even 

As for the even positions (rungs) the RR's ( ^9|) of the 
states are second order, we have to compute directly the 
values dii{R,k) and dii^i{R,k) {k = 1,2). Clearly we 
also have, 



dR-i{R,k) = 



(78) 



Now, as R is even, let us set R = 2r. Then, using the 
RR's we find. 



d2r{2r, k) — Z^r-l + ^r^\}2r-\ fe^2r- 



(79) 



In the derivation of these RR's we need to compute the 
following reduced density matrix elements, 

= ^1r\d2r.k\^2r) = 1 



{{hn hr~l)\d2r,k\{l2r, hr-l)) 



(80) 

-'2r.2r-l.k — \ V'^r j i2r- i y |"2r,/c | Vt2r j i2r- i ;/ (81) 

Using the expressions of the fluctuation states 
lihr, hr-i)) we find that the only nonvanishing elements 
are as follows, 

srCf^) _ 1 xi'^) — 1 A^'^) 



2;fc ~ 2 ' "l,4:fe ~ 2 ' "l,6;fc 



"2,l;k — 2' "44:fc ~ 2' "6,l;fc 



-'3,4;*: 



1, (5, 



(e) 



(82) 
(83) 
(84) 
(85) 



where we are following the notation for the site/rung 
states, {1,3,5} - {|o>,| i), | T>}, {2,4,6} = 
{|5),|Tt),|T,)}. 

Likewise we proceed with the other initial value 
(ii?,+i(i?) and we find. 



d2r+i(2r, k) = d2r{2r, k) + 

«r''2r+l,2r.fc^2r-l + Or''2r+l,2r,2r-l,fc^2r-2 



(86) 



In the derivation of these RR's we need to compute the 
following reduced density matrix elements, 

'^2r+l,2r,2r-l,fe = ((^2r+l , ^2r , ^2r-l ) M2r,fc | (^2r+l , ^2r , ^2r- 1 )) 

(87) 

Using the expressions of the fluctuation states 
|(Z2r+i, hr, hr-i)) wc find that the only nonvanishing el- 
ements are as follows. 



-'3,4,3;fc 



1, (5, 



5,6,5;fe 



APPENDIX B: SPECTRUM OF THE t-j MODEL 
ON THE PLAQUETTE 

The t — J Hamiltonian on the plaquette has the sym- 
metry group -D4 of the square. This implies that the 
eigenstates can be classified with the irreps of D4, to- 
gether with the number of holes h and the total spin 
S. D4 has 5 irreps, four of which are one-dimensional 
Ai, Bi, A2, B2 and one two-dimensional irrep E. From 
the character table of this group one sees that the irreps 
Ai, Bi have, in our terminology, even parity in both di- 
agonals, while the irreps A2 , B2 are odd. The parity of 
the states in the irrep E is (+, — ) and (— , -I-) for the two 
diagonals. In Table 3 we show the analytic expression of 
the energy and the value for the cases J = and J = 0.5 
(t = 1). In Fig. 33 we show a plot with all the energies 
as functions of J/t. Observe the crossover between the 
lowest energy states at J/t ^ 1.37. 
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Table 4. Here h is the number of holes, S the total spin, 
The third column denotes the irrep of the D4 group. We 
give the values of the energy for J = and 0.5 ( t = 1). 



1 



"3,4,l;fe ~ 2 



- 1 

' "5,6,l;fe ~ 2 



"l,4,3;fe 



1 S^'^ 
2' "l,5,6;fc 



(e) _ 1 
2 



(89) 
(90) 
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FIG. 33. Plot showing the energy levels for a tj model on 
the 2x2 plaquette as a function of the coupling J/t. The 
energy levels are classified according the D4 symmetry group 
of the square as shown in table 4. 



This suggests that there must be a direct way to re- 
late the valence bond construction of the spin 1 AKLT 
stat e |T^] and the dimer-RVB picture of the 2-leg lad- 
der | |28| , [l6| . We shall show that this is indeed possible 
through the plaquette construction of the 2-leg ladder, 
which is shown diagrammatically in Figs. 34 and 35. 

In Fig. 34 the 2-leg ladder is split into plaquettes con- 
nected by two links. We have generalized somewhat the 
medial construction in this case, since two plaquettes are 
allowed to have more than one common link. The inter- 
esting point is that the filling factor x of the ladder and 
the one of the extended model, related in exactly 

the same manner as the 2D lattices (see Eq.(Q)). Thus 
X = corresponds to Xp = 1/2. Let us assume for the 
moment that each plaquette has spin 1. Then the effec- 
tive interaction describing the coupling of these spins will 
be an AFH model. This suggests that we can construct a 
valence bond state to approximate the plaquette ground 
state by drawing nearest neighbors bonds among the ele- 
mentary spins between plaquettes, as shown in the upper 
part of Fig. 35. Now, if we project out of this state any 
5 = components of each plaquette, we get the AKLT 
state, where each plaquette is a pure spin S = I. If, 
instead, we Gutzwiller project the td links, as shown in 
Fig. 35, we get the dimer-RVB state proposed in refer- 
ences MM. 



APPENDIX C: PLAQUETTE DERIVATION OF 

THE EQUIVALENCE BETWEEN THE 
HALDANE STATE AND THE RVB STATE OF 
THE 2-LEG SPIN LADDER 

There has been some debate in the past as to whether 
the Haldane state of the spin 1 chain is in the same phase 
as the GS of the AFH 2-leg ladder. The general consensus 
is that they both belong to the same universality class 
characterized by a spin gap, finite spin correlation length 
and non vanishing string order parameter. The DMRG 
study of Ref. |p7| demonstrated a continuous mapping 
between these systems, and pointed out the equivalence 
between a dimer-RVB state on a composite spin model 
(which is a ladder model with some extra hopping terms) 
and the AKLT model. 



O •k^O o 

• o > o 



o 








f 








FIG. 35. (Below) A dimer-RVB state on a rectangular 2-leg 
ladder obtained as the projection (above) of a doped Xp = 1/2 
valence bond state on the decorated (dual) lattice associated 
to the rectangular ladder. 



FIG. 34. Plaquette construction (small interior squares 
plus dashed lines) of the rectangular ladder with n; = 2 legs. 



Hence, the plaquette model acts as intermediate sys- 
tem, for which different projections generate either the 
AKLT state of the spin-1 chain or the dimer-RVB state 
of the two-leg ladder. 
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